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■ Abstract. I describe very briefly the new libraries of empirical spectra of 

stars covering wide ranges of values of the atmospheric parameters T e //, log 
C^l , g, [Fe/H], as well as spectral type, that have become available in the recent 

£ . past, among them the HNGSL, MILES, UVES-POP, ELODIE, and the IndoUS 

libraries. I show the results of using the IndoUS and the HNGSL libraries, as 
well as an atlas of theoretical model atmospheres, to build population synthesis 
models. These libraries are complementary in spectral resolution and wavelength 
fT) ' coverage, and will prove extremely useful to describe spectral features expected 

in galaxy spectra from the NUV to the NIR. The fits to observed galaxy spectra 
using simple and composite stellar population models are discussed. 
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1. Introduction 

o . 

Bruzual fc CharlotJ ((20031 ) . hereafter BC03, have examined in detail the advan- 
tages of using intermediate resolution stellar spectra in population synthesis and 
' galaxy e volution models. In th ese models BC03 use the STELIB library com- 

piled by iLe Borgne et al. ( 20031 ). This library contains observed spectra of 249 



S3 



stars in a wide range of metallicities in the wavelength range from 3200 A to 9500 



A at a resolution of 3 A FWHM (corresponding to a median resolving power 
of A/ A A ~ 2000), with a sampling interval of 1 A and a signal-to-noise ratio of 
typically 50 per pixel. The BC03 models reproduce in deta il typical galaxy 



spectr a extracted from the SDSS Early Data Release, e.g., IStoughton et al 



X ' (|2002l ). From these spectral fits one can constrain physical parameters such 



as the star fo r matio n history, metallicity, a nd du s t content of galaxies, e.g. 
Heavens et~al] (|2004l ); ICid-Fernandes et al.l (|2005l ); iMateu. Magris. fc Bruzual 
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(2005). The medium resolution BC03 models also enable accurate studies of 



absorption-line strengths in galaxies containing stars over the full range of ages 
and can reproduce simultaneously the observed strengths of those Lick indices 
that do not depend strongly on elemental abundance ratios, provided that the 
observed velocity dispersion of the galaxies is accounted for properly, and offer 
the possibility to explore new indices over the full optical range of the STELIB 
atlas, i.e. 3200 A to 9500 A. To extend the spectral coverage in the models 
bey ond the s e limi ts, BC03 recur to other libraries. For solar metallicity models, 
the lPicklesI (|l998l ) library can be used to extend the STELIB spectral coverage 



down to 1150 A in the UV end and up to 2.5 /jm at the red end, with a sampling 
interval of 5 A pixel -1 and a median resolving power of (A/ A A ~ 500). The UV 
spectra in the iPicklei ( 19981 ) atlas are based on IUE spectra of bright stars. At 



all metallicities the BC03 models are extended down to 91 A in the UV side and 
up to 160 /jm at the other end using the theoretical model atmospheres included 
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in th e BaS eL series of libraries c ompiled bv lLejeune. Cuisinier. Sz Buserl ( 19971 . 
1998), and IWestera et al.1 ((2002;), but at a resolving power considerably lower 



than for STELIB (A/ A A 200 - 500). 

Given the success in reproducing observed galaxy spectra with the BC03 
synthesis models, it is important to explore models built with libraries of higher 
spectral resolution and which improve upon STELIB on the coverage of the 
HRD by including at all metallicities a broader and more complete distribution 
of spectral types. This goal is now possible thanks to several compilations of 
stellar spectra that have become available in the last few years. In this paper I 
describe briefly the properties of these libraries and show applications of galaxy 
models built using some of them. The full implementation of the new libraries 
in the population synthesis models is in preparation by Bruzual & Chariot. 



2. Stellar Libraries 



A relatively large large number (5) of libraries containing medium to high spec- 
tral resolution observed spectra of excellent quality for large numbers of stars 
(hundreds to thousand) have become available during the recent past. One of 
the main objectives of the observers who invested large amounts of time and 
effort assembling these data sets was to build libraries suitable for population 
synthesis. In this respect, the stars in the libraries have been selected to pro- 
vide broad coverage of the atmospheric parameters T e jj, log g, [Fe/H], as well 
as spectral type, throughout the HRD. In parallel to this observational effort, 
considerable progress has been made in the computation of theoretical model 
atmospheres at high spectral resolution for stars whose physical parameters are 
of interest for galaxy modeling. Below I describe very briefly the characteristics 
of each of these libraries which are relevant for population synthesis models. 



2.1. HNGSL 

The Hubble's New Generation Spectral Library ( Heap Sz Land 2003) contains 
spectra for a few hundred stars whose fundamental parameters, including chem- 
ical abundance, are well known from careful analysis of the visual spectrum. 
The spectra cover fully the wavelength range from 1700 A to 10,200 A. The 
advantage of this library over the ones listed below is the excellent coverage of 
the near-UV and the range from 9000 A to 10,200 A, which is generally noisy 
or absent in the other data sets. 



2.2. MILES 



The M edium resolution INT Library of Empirical Spectra (jSanchez-Blazquez et al 



2005), contains carefully calibrated and homogeneous quality spectra for 1003 



stars in the wavelength range 3500 A to 7500 A with 2 A spectral resolution and 
dispersion 0.9 A pixel -1 . The stars included in this library were chosen aiming 
at sampling stellar atmospheric parameters as completely as possible. 

2.3. UVES POP library 

The UVES Paranal Observatory Project (B agnulo et al. 20041 ). has produced a 



library of high resolution (R = A/AA ~ 80,000) and high signal-to-noise ratio 
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IndoUS vs. STELIB libraries: Metallicity distribution 
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Figure 1. Histogram showing the number of stars per 0.05 size bin in [Fe/H] 
for the IndoUS (dark line) and STELIB (hatched histogram) libraries. The 
values of [Fe/H] corresponding to the Padova 1994 and Padova 2000 sets of 
evolutionary tracks and the BaSeL 3.1 spectral atlas are shown as vertical 
lines at the top of the figure. 
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spectra for over 400 stars distributed throughout the HRD. For most of the 
spectra, the typical final SNR obtained in the V band is between 300 and 500. 
The UVES POP library is the richest available database of observed optical 
spectral lines. 

2.4. IndoUS library 

The IndoUS library ( Valdes et al .1120041 ) contains complete spectra over the entire 



3460 A to 9464 A wavelength region for 885 stars obtained with the 0.9m Coude 
Feed telescope at KPNO. The spectral resolution is ~ 1 A and the dispersion 
0.44 A pixel -1 . The library includes data for an additional 388 stars, but only 
with partial spectral coverage. 

2.5. ELODIE 

The ELODIE library is a stellar database of 1959 spectra for 1503 stars, ob- 
served with the echelle spectrograph ELODIE on the 193 cm telescope at the 
Observatoire de Haute Provence. The resolutio n of this library is R = 42, 00 
in the wavelength range from 4000 A to 6800 A dPrugniel & Soubiranll2001allbh . 



This library has been updated, extended, and used by iLe Borgne et al.l ( 20041 ) 
in the version 2 of the population synthesis code PEGASE. 

2.6. High-spectral resolution theoretical libraries 

There are several on-going efforts to improve the existing grids of theoreti- 
cal model atmospheres including the computation of high resolution theoret- 
ical spectra for stars whose phys ical paramete r s are o f interest for popula - 
tion synthesis. See, for e xample, Coelho et al.l (|2005l ); iBertone et al.l ( 20041 ) ; 
Murphv k Meiksinl (|2003 l : iPeterson et al.1 (j2004j ). as well as the papers by Chavez, 
Bertone, Buzzoni & Rodriguez-Merino; Gonzalez-Delgado & Cerviho; and Mu- 
nari & Castelli in this conference. 



3. Using the New Libraries in Population Synthesis Models 

In BC03 the 'standard' reference model represents a simple st ellar population 
(SSP) computed using the Padova 1994 evolutionary tracks, the lChabrierj (2003) 



IMF truncated at 0.1 and 100 M , and either the STELIB, the Pickles, or the 
BaSeL 3.1 spectral libraries (see BC03 for details). Below I show the behavior 
of the standard reference model com puted with the IndoUS and the HNGSL 
libraries, and the lCoelho et all ( 20051 ) atlas of theoretical model atmospheres. 



3.1. IndoUS library 

Figure 1 shows a histogram of the number of stars per 0.05 size bin in [Fe/H] 
for both the IndoUS and the STELIB libraries. The gain in the number of stars 
available at most metallicities in the IndoUS library with respect to STELIB 
shows clearly in this plot. This should translate in a better sampling of most 
stellar types at the relevant positions in the HRD for the metallicities corre- 
sponding to the different evolutionary tracks, also shown in Figure 1. Models 
built using the IndoUS library are compared to the STELIB and BaSeL 3.1 
models in Figures 2 to 5. As expected, spectral features are more numerous and 
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Figure 2. BC03 standard SSP model spectra for solar metallicity at 12 Gyr 
in the wavelength range from 3500 A to 9500 A computed with the IndoUS 
(top line) and the STELIB (bottom line) libraries. The spectra have been 
scaled and shifted in the vertical direction for clarity. 
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Indo-US at 12 Gyr, Z = 0.004, 0.008, 0.02, 0.05 
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Figure 3. BC03 standard SSP model spectra for Z = 0.004, 0.008, 0.02, 
and 0.05 (top to bottom) at 12 Gyr in the wavelength range from 3500 to 9500 
A built using the IndoUS stellar library. The spectra have been normalized 
and shifted arbitrarily in the vertical direction for clarity. 
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STELIB at 12 Gyr, Z = 0.004, 0.008, 0.02, 0.05 
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Figure 4. BC03 standard SSP model spectra for Z = 0.004, 0.008, 0.02, 
and 0.05 (top to bottom) at 12 Gyr in the wavelength range from 3500 to 9500 
A built using the STELIB stellar library. The spectra have been normalized 
and shifted arbitrarily in the vertical direction for clarity. 
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BaSeL 3.1 at 12 Gyr, Z = 0.004, 0.008, 0.02, 0.05 
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Figure 5. BC03 standard SSP model spectra for Z = 0.004, 0.008, 0.02, 
and 0.05 (top to bottom) at 12 Gyr in the wavelength range from 3500 to 9500 
A built using the BaSeL 3.1 spectral atlas. The spectra have been normalized 
and shifted arbitrarily in the vertical direction for clarity. 
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more clearly seen in the IndoUS model than in the other two models. Colors 
measured from the continuum flux in these three spectra are all very similar. 

Figures 6 and 7 show the normalized residual distributions resulting from 
simultaneous fits of the strengths of several indices in the spectra of 2010 galaxies 
with S/N me( j > 30 in the 'main galaxy sample' of the SDSS EDR using BC03 
models built with the STELIB and the IndoUS libraries, respectively. The 
highlighted frames indicate the seven indices used to constrain the fits. The 
stellar velocity dispersion of the models are required to be within 15 km s _1 
of the observed ones. Each panel shows the distribution of the fitted index 
strength / minus the observed one I obs , divided by the associated error crj. For 
reference, a dotted line in each panel indicates a Gaussian distribution with unit 
standard deviation. The shaded histograms show the contributions to the total 
distributions by galaxies with ay > 180 km s , corresponding roughly to the 
median stellar velocity dispersion of the sample. Figure 6 is identical to Figure 
18 of BC03. It is apparent from a comparison of these two figures that the CNi, 
CN2, Ca4227, G4300, TiOi, and TiC"2 indices measured in this galaxy sample 
are reproduced more closely by the IndoUS models than the STELIB models. 
The Mgi, Mg2, Mgb, and NaD indices are equally off in both sets of models. 
The negative residuals indicate that the model values are below the observed 
ones, due to the lack of stars with enhanced a element abundance in these stellar 
libraries. The distributions for H/3 and H^ in the IndoUS models are shifted 
toward positive and negative residual values, respectively, with respect to the 
STELIB models. The four [MgFe] indices shown in the bottom row of Figures 
6 and 7 are all shifted toward negative residuals in the IndoUS models. 



Flux calibration problems in the IndoUS library. In Figure 8 I compare the 
spectra of 6 of the more than 500 stars in common between the IndoUS and the 
MILES libraries. These plots illustrate a flux calibration problem affecting some 
of the IndoUS spectra. Figure 9 shows spectra of 6 of the stars in common in the 
IndoUS, the MILES, and the STELIB libraries. For all the stars in common in 
the three libraries, the shape of the MILES and STELIB spectra always match, 
whereas if there is a discrepancy it can always be attributed to the IndoUS 
spectrum, as in the case of the star HD122563 in the upper right hand side panel 
of Figure 9. The IndoUS spectra of stars with flux calibration problems should 
not be used in population synthesis unless they can be re-calibrated correctly. 
In Figure 10 I indicate a procedure that uses the BaSeL 3.1 spectrum for the 
values of T e ff, log g, and [Fe/H] of the problem star to define the correct shape 
of the IndoUS spectrum. This procedure assumes that the physical parameters 
of the star are known accurately. 



3.2. HNGSL 



I hav e also computed the BC03 standard model using the HNGSL (|Heap fc Land 
20031 ) instead of STELIB to represent the stellar spectra. Figure 11 shows clearly 



the advantages of using the HNGSL to study the near UV below 3300 A. Spectral 
features, lines and discontinuities, are much better defined in the HNGSL spectra 
than in the IU E spectr a (used in the Pickles library) and the BaSeL 3.1 atlas 
(see BC03 and Bruzual ( 20041 ) for details). However, colors measured from the 
continuum flux in these three spectra are very similar. The higher spectral 
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Figure 6. Simultaneous fit of the strengths of several indices in the spectra 
of 2010 galaxies with S/N moc j > 30 in the 'main galaxy sample' of the SDSS 
EDR using BC03 models built with the STELIB library. The highlighted 
frames indicate the seven indices used to constrain the fits. The stellar ve- 
locity dispersion of the models are required to be within 15 km s _1 of the 
observed ones. Each panel shows the distribution of the fitted index strength 
jfit mmus t ne observed one 7 obs , divided by the associated error 07. For ref- 
erence, a dotted line in each panel indicates a Gaussian distribution with unit 
standard deviation. The shaded histograms show the contributions to the 
total distributions bv ealaxics with ov > 180 km s . corrcsDondine rouehlv 
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Figure 7. Simultaneous fit of the strengths of several indices in the spectra 
of 2010 galaxies with S/N moc j > 30 in the 'main galaxy sample' of the SDSS 
EDR using BC03 models built with the IndoUS library. The highlighted 
frames indicate the seven indices used to constrain the fits. The stellar ve- 
locity dispersion of the models are required to be within 15 km s _1 of the 
observed ones. Each panel shows the distribution of the fitted index strength 
jfit mmus t ne observed one 7 obs , divided by the associated error 07. For ref- 
erence, a dotted line in each panel indicates a Gaussian distribution with unit 
standard deviation. The shaded histograms show the contributions to the 
total distributions bv ealaxics with ov > 180 km s . corrcsDondine rouehlv 
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Figure 8. Comparison of the sed's of 6 stars present in the IndoUS and the 
MILES stellar libraries in the wavelength range AA 3500 - 7400 A in common 
in these libraries. In each panel the black line is the IndoUS sed, the light 
gray line is the MILES sed, and the dark gray line represents the difference 
(IndoUS - MILES). The sed's have been scaled arbitrarily in the flux scale. 
The star identification is indicated in the upper left corner of each panel. 
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Figure 9. Comparison of the sed's of 6 stars present in the IndoUS, the 
STELIB, and the MILES stellar libraries in the wavelength range A A 3500 
- 7400 A in common in these libraries. In each panel the black line is the 
IndoUS sed, the light gray line is the MILES sed, and the dark gray the 
STELIB sed. The sed's have been scaled arbitrarily in the flux scale. The 
star identification is indicated in the upper left corner of each panel. 



G. Bruzual 



HD231170, G2IV. [Fe/H] - -2.19, T rf[ - 4500, log g - 0.90 
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Figure 10. Spectral energy distributions of the star HD221170. (a). The 
MILES sed extends up to 7400 A and the IndoUS sed up to 9500 A. (b). 
Smoothed version of the IndoUS sed for this star and of the BaSeL 3.1 model 
for T e ff = 4500K, log g = 1, and [Fe/H] = -2, corresponding to the tabulated 
values for HD221170. The residual line in the bottom part of the panel is the 
difference between the original and the smoothed IndoUS sed's. (c). Result of 
adding the smoothed BaSeL sed to the difference of the original and smoothed 
IndoUS sed's shown in frame (b). 
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resolution of STELIB above 3300 A compared to HNGSL is clearly seen in this 
figure. 

The spectra in Figure 12 are useful to study the different behavior of models 
computed with different libraries in the region around and above 9000 A. There 
is a clear difference in the continuum level predicted by the HNGSL and STELIB 
models (lowest level) compared to the level predicted by the BaSeL 3.1 and the 
Pickles library models (highest level). The HNGSL model is to be preferred, 
given the higher SNR of this spectrum in this region. However, spectral features 
below 8500 A are more clearly seen in the higher resolution STELIB model. 



3.3. Theoretical model atmospheres 

Figures 13 and 14 show the BC03 standard SSP model sed at 12 Gyr computed 
with the high resolution theoretical model atmospheres from the IAG collabo- 
ration for [Fe/H] = 0, [ct/Fe] = 0, in the wavelength range from 0.3 to 18 (j,m 
(ICoelho et alJ 120051 ). . n Figure 14 one can clearly see those regions of the spec- 



trum where the theoretical atmosphere models produce different results from the 
empirical library models, e.g., above 9000 A. In Figure 15 I compare the STELIB 
and IndoUS models with the IAG collaboration model degraded in spectral res- 
olution to match the empirical spectra. This figure shows a remarkable degree 
of agreement between the empirical and theoretical models when the difference 
in resolution is duly taken into account. This figure replaces the one I showed 
in my oral presentation which contained an involuntary error in the procedure 
used to downgrade the spectral resolution of the theoretical model. 



4. Reproducing Galaxy SEDs 

4.1. SSP fits 

By means of a standard least-squares technique it is possible to select the age 
and stellar velocity dispersion at which a given SSP model reproduces most 
closely a given observed galaxy spectrum. Figures 16 to 18 show remarkably 
good fits to the continuum spectrum of three galaxies with different ages and 
velocity dispersions, and in different wavelength intervals. The BC03 standard 
SSP model computed with the IndoUS stellar library was used in these fits. The 
stellar velocity dispersion a indicated in each figure and applied to the model 
spectrum is also derived by the fitting algorithm. The residuals (observed - 
model) shown at the bottom of the plots are quite flat over the whole spectral 
range. 

4.2. Non-parametric CSP fits 

An alternative approach for studying the stellar populations present in a galaxy, 
based on the fact that galaxies are thought to be more closely described by 
a composite stellar population (CSP) rather than by an SSP, consists in de- 
riving the galaxy star formation history (SFH) by means of a non-parametric 
CSP fit to the galaxy sed. In a non-parametric fit the SFH is not assumed 
to be described by an analytic function, e.g., an exponentially decaying star 
formation rate with characteristic e— folding time r. Instead, all the spectra 
defining the evolution of an SSP of one or more metallicity values are allowed 
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1 Gyr, Z G , Chabrier IMF, Padova 1994 
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Figure 11. BC03 standard SSP model spectra for solar metallicity at 1 Gyr 
in the wavelength range from 1600 A to 3800 A. The spectrum on top (thick 
black line) is built using the HNGSL, the one in the middle (gray line) uses 
the Pickles library and STELIB, and the bottom one (thin black line) the 
BaSel 3.1 library. The spectra have been normalized at 2000 A and shifted 
arbitrarily in the vertical direction for clarity. 



Studying stellar populations at high spectral resolution 



17 



12 Gyr, Z , Chabrier IMF, Padova 1994 




4000 6000 8000 10 4 



wavelength (A) 

Figure 12. BC03 standard SSP model spectra for solar metallicity at 12 
Gyr in the wavelength range from 3500 A to 10,500 A. The black solid line 
(reaching up to 10,200 A) represents the spectrum built using the HNGSL, 
the dark gray solid line (noisy above 9000 A) uses STELIB, the light gray solid 
line uses STELIB up to 9000 A and the Pickles library at longer wavelengths, 
the black dashed line uses the BaSel 3.1 library. The spectra have been 
normalized at 5000 A. 
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Figure 13. BC03 standard SSP model sed at 12 Gyr computed with the 
high resolution theoretic al model atmosphe res from the IAG collaboration for 
[Fe/H] = 0, [a/Fe] = (jCoelho et al.ll2005l) in the wavelength range from 0.3 
to 18 fim. 
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Figure 14. BC03 standard SSP model sed at 12 Gyr computed with the 
high resolution theoreti cal model atmosphe res from the IAG collaboration for 
[Fe/H] = 0, [ot/Fe] = (jCoelho et alJ2005f ) in the wavelength range from 3500 
to 9500 A (gray line) downgraded to the spectral resolution of the IndoUS 
library. For comparison, the corresponding model computed with the IndoUS 
stellar library is shown as a black line. 
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Figure 15. BC03 standard SSP model sed at 12 Gyr computed with the 
high resolution theoretical model atmosp heres from the IAG collaboration for 
[Fe/H] = 0, [a/Fe] = 0.0 (jCoelho et alj|2005l ). in the wavelength range from 
5000 to 5200 A containing the Mgi, Mg 2 , and Mgb features. The thin lines 
show the IAG model downgraded to the spectral resolution of the STELIB 
and IndoUS libraries, shown as heavy lines. The spectra have been scaled and 
shifted arbitrarily in the vertical direction. 
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to be selected individually by a spectral fitting algorithm. The SFH for the 
galaxy is re built from the known age and m etallicity of each of the selected 
spectra. See lMateu. Magris. &: Bruzuall ( 20051 ) for a description of GASPEX, a 
non-parametric CSP fitting procedure based on t he Non-Negative Least Squares 
(NNLS) algorithm of lLawson and Hanson! (jl974l ). 

The SFH recovered by this type of algorithm is highly dependent on the 
signal-to-noise ratio and the wavelength range covered by the problem spectrum. 
Figures 19 and 20 show the results of applying the GASPEX algorithm to a 
problem spectrum of known SFH. The problem spectrum corresponds to an old 
stellar population of age 6 Gyr in which a second burst of star formation of i 
the intensity of the first burst occurs 360 Myr ago. In these figures, the two 
vertical lines at the corresponding ages represent schematically this SFH. The 
histograms with error bars in Figure 19 represent the SFH recovered by the 
GASPEX algorithm when noise in varying amounts is added to the problem 
spectrum. The value of the signal-to-noise ratio (SNR) used is indicated in each 
frame. Values of SNR above 10 are required to recover the input SFH. Similarly, 
Figure 20 shows that the broader the wavelength range covered by the problem 
spectrum, the more likely it is that we recover the input SFH. For the example 
in Figure 20, the problem spectrum is the same one used in Figure 19 for a 
signal-to-noise ratio of 20. The full coverage of the wavelength range from the 
UV to the IR is required to recover the known SFH. 

Figure 21 shows the SFH and cumulative mass vs. age function recovered 
by the GASPEX algorithm for a selection of galaxies from the Near Field Galaxy 
Survey ( Jansen et al. 2000] ) . The thin lines represent galaxies with significant 



recent star formation, as indicated by the emission lines in their spectra. It is 
apparent from this plot that these galaxies assemble their stars at a later time 
than the early- type galaxies represented by the heavy lines in Figure 21. 



4.3. SSP vs. CSP fitting 

Figure 22 shows the best SSP model fit to the spectrum of an early type galaxy 
resulting from co-adding the spectra of several galaxies from the SDSS until 
reaching a signal-to-noise ratio of 250. This spectrum, kindly provided by Jarle 
Brinchmann, is shown as a gray line in the wavelength range from 3600 to 
7500 A. The BC03 standard SSP model computed with the STELIB stellar 
library for solar metallicity at an age of 9.25 Gyr is shown in black. The stellar 
velocity dispersion a = 280 km s _1 applied to the model spectrum is also derived 
by the fitting algorithm. At the bottom of the figure the line marked SSP 
represents the residuals (observed -model) for the SSP model fit shown in the 
figure. The line marked CSP represents the residuals obtained when the fit is 
performed by the non-parametric GASPEX CSP fitting algorithm described by 
Mateu. Magris. Sz Bruzual (2005). The line marked S-C shows the difference 



between the SSP and CSP solutions (shifted down arbitrarily for clarity). 

The top frame of Figure 23 shows the cumulative squared residuals (oc x 2 ) 
for the SSP and CSP fits shown in Figure 22. The CSP fit reaches a lower x 2 
value. The improved fit over the SSP fit is achieved by the GASPEX algorithm 
by assuming that some star formation happened in this galaxy at ages younger 
than the SSP age of 9 Gyr. The vertical lines show the position of the central 
bands defining the CN1-CN2, Mgi-Mg2-Mgb, and NaD Lick indices, and the 
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Ha line. The lower residuals in the CSP fit are possible because the younger 
population fills in the deficiency of the SSP model at the bands containing the 
Mg and Na features, i.e. the steep increase in \ 2 seen at these wavelengths in 
the SSP fit are not present in the CSP fit. On the other hand, the CN features 
are reproduced best by the SSP model. The Ha emission line is not included in 
neither of the models. The same is true for the apparent features in the observed 
spectrum at A > 6700 A. Beyond Ha the two x 2 functions run almost parallel. 

The middle frame of figure 23 shows the percentage contribution of the old 
and a very young population (100 Myr) to the total spectrum of this galaxy in 
the CSP solution as a function of wavelength. The contribution of populations 
of other ages included in the GASPEX solution is, added all together, below the 
5% level at A < 4500 A, and is even less at longer wavelengths. The spectra 
corresponding to the old, young, and the rest of the stellar populations in the 
GASPEX solution are shown in the bottom frame of this figure. The latter 
contribute essentially zero flux to the total sed. The heavy line in this frame 
represents the CSP solution, i.e., the addition of the individual spectra. 

The fact that the best GASPEX solution has a lower x 2 an d requires recent 
star formation does not necessarily imply that it represents a more realistic 
solution than the SSP fit from the astrophysical point of view. However, the 
GASPEX non-parametric CSP solution for this galaxy es reminiscent of the 
E+A phenomenon, quite common in E galaxies. The fact that the CSP solution 
seems able to detect the E+A phenomenon makes this technique appealing. 

Figure 24 shows a detailed comparison of the fit to one of J. Brinchmann 
co-added spectra in specific wavelength regions. The residuals in the second, 
fourth and sixth panels are expressed in units of the standard deviation a of the 
co-added fluxes. The signal-to-noise ratio in these spectra is so high, and a so 
low, that the very small differences between the problem and solution spectra 
translate into very large residuals when measured in units of a. Thus, very 
careful analysis of the fits should be performed when deriving galaxy SFHs. 



5. Conclusions 

New libraries of empirical spectra of stars covering wide ranges of values of 
the atmospheric parameters T e jj, log g, [Fe/H], as well as spectral type, that 
have become available recently, are complementary in spectral resolution and 
wavelength coverage. 

Models built using the IndoUS and the HNGSL empirical stellar libraries, 
and the IAG collaboration theoretical model atmospheres, show that this sort 
of library will prove extremely useful in describing spectral features expected in 
galaxy spectra of various ages and metallicities from the NUV to the NIR. 

Models build with the IndoUS library increase the number of Lick indices 
that are correctly reproduced by the models. 

The use of population synthesis models to derive galaxy star formation 
histories by fitting galaxy spectra have become very popular in recent years. The 
reader should keep in mind that different solutions are provided by the simple 
minimum \ 2 SSP fit and the more complex non-parametric CSP fit. Even if 
the CSP fits result in a lower value of x 2 > the solutions are not necessarily more 
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sound from astrophysical grounds. The fact that the CSP solution seems able 
to detect the E+A phenomenon makes this technique appealing. 

Complete sets of models that use the libraries mentioned in §2 are been 
built and will be discussed in a coming paper by Bruzual & Chariot. 

Acknowledgments. I thank Jarle Brinchmann for providing me with his 
sample of co-added SDSS spectra. 
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Figure 16. Best SSP model fit to the spectrum of the elliptical galaxy NGC 
221 in the wavelength range from 3600 to 5600 A. The BC03 standard SSP 
model computed with the IndoUS stellar library for solar metallicity at an age 
of 3.5 Gyr is shown in black. The stellar velocity dispersion a = 150 km s -1 
applied to the model spectrum is also derived by the fitting algorithm. The 
observed spectrum is shown in gray and was kindly provided by P. Sanchez- 
Blazquez. The residuals (observed - model) are shown at the bottom of the 
plot in the same vertical scale. 
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Figure 17. Best SSP model fit to the spectrum of the early-type galaxy 
SDSS-267-110 in the wavelength range from 3750 to 8800 A. The BC03 stan- 
dard SSP model computed with the IndoUS stellar library for solar metallicity 
at an age of 5 Gyr is shown in black. The stellar velocity dispersion a = 150 
km s _1 applied to the model spectrum is also derived by the fitting algorithm. 
The SDSS-EDR observed spectrum is shown in gray. The residuals (observed 
- model) are shown at the bottom of the plot in the same vertical scale. 
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Figure 18. Best SSP model fit to the continuum spectrum of the star 
forming galaxy SDSS-52365-0879-580 in the wavelength range from 3500 to 
9500 A. The BC03 standard SSP model computed with the IndoUS stellar 
library for solar metallicity at an age of 6 Gyr is shown in black. The stellar 
velocity dispersion a = 100 km s _1 applied to the model spectrum is also 
derived by the fitting algorithm. The SDSS-EDR observed spectrum is shown 
in gray. The wavelength ranges containing emission lines in the observed 
spectrum are not included in the fit. 
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Figure 19. Recovery of a known star formation history as a function of 
the signal-to-noise ratio of the problem spectrum. The problem spectrum 
corresponds to an old stellar population of age 6 Gyr in which a second burst 
of star formation of ^ the intensity of the first burst occurs 360 Myr ago. 
The two vertical lines at the corresponding ages represent schematically this 
star formation history. The histograms with error bars represent the star 
formation history rec overed by the GASPEX non-param etric CSP fitting al- 
gorithm described by iMateu. Magris. fe; Bruzual (|2005f) when noise is added 
to the problem spectrum. The value of the signal-to-noise ratio (SNR) used 
is indicated in each frame. 



28 



G. Bruzual 



I From 5000 A to 7000 A 



0.6 
0.4 




0.1 



10 0.1 

Age / Gyrs 



10 



Figure 20. Recovery of a known star formation history as a function of 
wavelength coverage. The problem spectrum is the same one used in Figure 
19 for a signal-to-noise ratio of 20. The histograms with error bars represent 
the star formation history re covered by the GASPEX non-pa rametric CSP 
fitting algorithm described bv lMateu. Magris. fc Bruz ual ( 2005) as a function 
of the wavelength range used in the fit, as indicated in each frame. 



Studying stellar populations at high spectral resolution 



29 




Figure 21. Star formation history and cumulative mass vs. age recovered 
by the GASPEX non-parametric CSP fitting algorithm described by Mateu, 
Magris, & Bruzual (2005) for a selection of galaxies from the Near Field 
Galaxy Survey (Jansen et al. 2000). The gray lines represent galaxies with 
significant recent star formation, as indicated by the emission lines in their 
spectra. 
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Figure 22. Best model fit to the spectrum of an early type galaxy resulting 
from co-adding the spectra of several galaxies from the SDSS until reaching a 
signal-to-noise ratio of 250. This spectrum, kindly provided by Jarlc Brinch- 
mann, is shown as a light line in the wavelength range from 3600 to 7500 A. 
The BC03 standard SSP model computed with the STELIB stellar library for 
solar metallicity at an age of 9.25 Gyr is shown as a heavy line. The stellar 
velocity dispersion a — 280 km s _1 applied to the model spectrum is also 
derived by the fitting algorithm. At the bottom of the figure the line marked 
SSP represents the residuals (observed -model) for the SSP model shown in 
the figure. The line marked CSP represents the residuals obtained when the 
fit is performed by the GASPEX non-parametric CSP fitting algorithm de- 
scribed by Mateu, Magris, & Bruzual (2005). The line marked S-C shows the 
difference between the SSP and CSP solutions. The last two lines have been 
shifted down for clarity. 
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Figure 23. Top frame. Cumulative squared residuals for the SSP and CSP 
fits shown in Figure 22. The vertical lines show the position of the central 
bands defining the CNi-CN 2 , Mgi-Mg 2 -Mgb, and NaD Lick indices, and the 
Ha line. Middle frame. Percentage contribution of the old and a very young 
population (100 Myr) to the total spectrum of this galaxy in the CSP solution 
as a function of wavelength. The contribution of populations of other ages 
included in the GASPEX solution is, added all together, below the 5% level 
at A < 4500 A, and is even less at longer wavelengths. Bottom frame. Spectra 
corresponding to the old, young, and the rest of the stellar populations in the 
GASPEX solution. The latter contribute essentially zero flux to the total sed, 
represented by the heavy line. 
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Figure 24. Detailed comparison of the fit to one of J. Brinchmann co-added 
spectrum in specific wavelength regions. The residuals in the second, fourth 
and sixth panels are expressed in units of the standard deviation a of the 
co-added fluxes. 



